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ABSTRACT 

We present the results of a high-resolution imaging survey of 12 brown dwarfs and very low mass 
stars in the closest (~145 pc) young (~5 Myr) OB association, Upper Scorpius. We obtained images 
with the Advanced Camera for Surveys/High Resolution Camera on HST through the f (U), 
F775W {i'), and F850LP {z') filters. This survey discovered three new binary systems, including one 
marginally resolved pair with a projected separation of only 4.9 AU, resulting in an observed binary 
fraction of 25 ± 14% at separations >4 AU. After correcting for detection biases assuming a uniform 
distribution of mass ratios for rris/mp > 0.6, the estimated binary fraction is 33 ± 17%. The binary 
fraction is consistent with that inferred for higher-mass stars in Upper Sco, but the separation and 
mass ratio distributions appear to be different. All three low- mass binary systems in Upper Sco are 
tight (< 18 AU) and of similar mass [nis/mp > 0.6), consistent with expectations based on previous 
multiplicity studies of brown dwarfs and very low mass stars in the field and in open clusters. The 
implication is that the distinct separation and mass ratio distributions of low-mass systems are set 
in the formation process or at very young ages, rather than by dynamical disruption of wide systems 
at ages >5 Myr. Finally, we combine the survey detection limits with the models of Burrows et al. 
(1997) to show that there are no planets or very low-mass brown dwarfs with masses > lOAfj at 
projected separations > 20 AU, or masses > 5Mj at projected separations > 40 AU orbiting any of 
the low-mass (0.04-0.10 Mq) objects in our sample. 

Subject headings: stars:binaries:visual — stars:low-mass, brown dwarfs — stars:pre-main sequence 



1. INTRODUCTION 

Brown dwarfs (BDs) are objects with masses inter- 
mediate between those of stars and planets. As such, 
studies of BDs potentially can offer a unique link be- 
tween theories of star formation and planet forma- 
tion. Several BD formation scenarios have been pro- 
posed, but current observational results indicate that the 
properties of BDs, includi ng their spatial distribution 
in star- forming regions ( Luhman 2004b; Briceno et al. 
[2002,.). velocity dispe rsion (jJocrgons .& Gucnthcr 2001: 
iWhite fc Bggj|_| 2003t). accr e tion and inner disk fre- 
quency llMuench et alJl200H iJavawardhan a et al.l l2003t 
iLiu et alJl2003t iMohanty et alJl 2005D. and evidence for 
an embedded phase (jWhite Hillenbrandl 12004"), are 
consistent with those of low-mass stars. Although this 
strongly supports the idea that BDs form in a fashion 
similar to stars, via the dynamical collapse and fragmen- 
tation of a cloud core, there is emerging evidence that 
the binary properties of brown dwarfs are fundamentally 
different from those of low-mass stars. 

Multiplicity sur veys of field T dwarfs ( Buraasscr e^^lj 
200311. L dwarfs JKoerner et aplTggft .T^eid et a,ljT20mr 
Close et all 1200^ IBouv et al.l 12003b IGizis et alJ Mm . 
late M dwarfs (iSiegl er et al.l 12005) and intermediate- 
age BDs in the Pleiades l|Martin et al. 2000) have 
found lower binary frequencies and smaller binary 
separations than for field G and early M dw arfs 
l|Duauennov fc MavoJlTOQlt iFischer fc Marcvlll992D . It 
is unclear if these binary properties reflect differences 
in the formation mechanism or the dynamical disrup- 
tion of wide substellar binaries. One way to distin- 
guish between these scenarios is by studying the prop- 



erties of young BD binaries, possibly before any signif- 
icant dynamical evolution has occurred. In this paper 
we present the results of an imaging multiplicity survey 
of br own dwarfs in the nearest (~145 pc) young (~5 
Mvr) l|Preibisch fc Zinneckeilll999D OB Association, Up- 
per Scorpius. 

2. OBSERVATIONS AND DATA REDUCTION 

Our targets were selected from the survey for low-mass 
members of Upper Sco by Ardila ct al. (200Q.). They 
identified 15 candidate members with spectral types of 
M5.5 or cooler, based on either low resolution spectra 
(10 objects) 01 R — Ic colors (5 objects). Three of these 
15 (USco 85, 114, 121) were determined to be likely non- 
members, based on low lithium abundance and radial 
veloc ities that are inconsistent with higher rn ass mem- 
bers HMuzerolle et 31.11200^ IWhite et al.il200l. The re- 
maining 12 were observed in this program. Subsequently, 
membership for all but the faintest member of this sam- 
ple, USco 137, has been confirmed with additional high 
resolution spectra l|Mohantv et al.ll2'005D . USco 137 has 
not yet been observed spectroscopically. Since our tar- 
gets include both very low mass stars and brown dwarfs, 
we hereafter refer to them as very low mass objects, or 
VLMOs. 

Our images were obtained with the Advanced Cam- 
era for Surveys/High Resolution Camera on the Hubble 
Space Telescope, which has a field of view of 26x29"aiid 
distortion-corrected pixel size of 27 mas pix^^. Observa- 
tions were made between July and September 2003 with 
the filters F555W (V), F775W (i'), and F850LP (z') at 
two dither positions near the center of the detector and 
with two exposures per position. Total integration times 
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were 510, 300, and 200 seconds, respectively. The F555W 
exposure times for the brightest objects {USco-55, 66, 
67, and 75) was reduced to 350 seconds to allow for ad- 
ditional short exposures in F775W and F850LP, which 
were close to the saturation limit in the full-length expo- 
sures. We chose the V band to maximize angular resolu- 
tion (diffraction limit 9resy = 58 mas) and the i' and z' 
bands to maximize sensitivity to very low mass compan- 
ions. 

The raw images were calibrated and distortion- 
corrected by the CALAC S pipeline during on-the-fly- 
reprocessing (IMack et al.l 120 03). Some cosmic rays re- 
mained, but their morphologies were substantially dif- 
ferent from stellar PSFs, so they were easily identified 
by visual inspection. 

Potential point sources were identified with the IRAF 
task DAOPHOT/DAOFIND, which found all local 
brightness maxima with a significance of > 5(t and a 
full-width at half-maximum (FWHM) near the expected 
value for the filter. We then measured aperture photom- 
etry and point-spread function (PSF) fitting photome- 
try for all objects in each field, and we report aperture 
photometry for all isolated objects and PSF photometry 
for all close binaries. PSF magnitudes were corrected to 
match aperture magnitudes based on results for the nine 
well-sampled objects that appeared isolated under visual 
inspection. Aperture photometry was carried out with 
the IRAF task DAOPHOT/PHOT with a 5-pixel aper- 
ture for faint objects and a 10-pixel aperture for bright 
objects, and our sky annulus had an inner radius of 200 
pixels and width of 15 pixels. We used the finite aperture 
corrections of iSirianni et alJ l)2005() . Point-spread func- 
tion (PSF) fitting photometry was carried out wi th the 
IRAF task DAOPHOT/ALLSTAR l|Stetsonlll987|) . The 
PSF for each filter was constructed from the 9 isolated 
VLMOs. Since all of the targets were located near the 
center of the chip and have similar temperatures and ex- 
tinctions, location and color effects should not be impor- 
tant. There was some variation in the PSF FWHM from 
target to target (±5%), which we attribute to small orbit- 
to-orbit changes in focus. However, the radial shape of 
the PSF was very similar for all objects except USco-109, 
which we will discuss in more detail in Section 3.1. 

Transformations to ground-based magnitudes {V , 
SDSS i', and SDSS z') were calculated with the IRAF 
task SYNPHOT/CALCPHOT, which convolves an in- 
put spectrum with transmission curves for HST's optics 
and filters or standard ground-based filters. Since SYN- 
PHOT does not include transmission curves for the i' and 
z' filters, which are defined at the United States Naval 
Observatory 40- in telescope, they were obtained from the 
website for the Sloan Digital Sky Survey^ and represent 
the filters and opti cs at the USN O-40 observing at 1.3 
airmasses iFukuaa wa et al.lll996f) . Based on transforma- 
tions determined for a set of M0-M8 dwarfs and M5- 
M8 giants from the Bru zual-Persson-Gunn -Stryker Spec- 
trophotometry Atlas HBruzual et al.l ll996). we find con- 
stant corrections that do not depend significantly on tem- 
perature or surface gravity: msss — my — —0.16 ± 0.03, 
mrr^-rm, = -h0.07±0.03, andmsso-m^' = -h0.03±0.03. 
The uncertainties are estimated from the standard devia- 
tion between all tested objects. The transformed magni- 

http: / / www.sdss.org/drl/ algorithms / standardstars/Filters / res 



tudes^ are listed in Table 1. The statistical uncertainties 
correspond to either the photon noise (for aperture pho- 
tometry; typically <0.01 magnitudes) or the goodness of 
fit (for PSF-fitting photometry; 0.02-0.09 magnitudes). 
Systematic uncertainties in the magnitude transforma- 
tions and aperture corrections are ~0.03 magnitudes. 

3. RESULTS 
3.1. New VLMO Binaries 

In Figure 1, we present contour plots of three candidate 
binaries (USco-55, 66, and 109) and the apparently single 
star USco-67 in the F555W, F775W, and F850LP filters. 
The USco-55 and USco-66 systems are clearly resolved. 
The USco-109 system is not obviously resolved and was 
not initially reported as a double-source by ALLSTAR, 
but the PSF appears to be marginally elongated in the 
-l-x direction relative to the single VLMO USco-67. 

One limitation in the ALLSTAR-based data reduction 
method is that binaries with very close (^A/Z?) separa- 
tions are often not identified, even when their combined 
PSF is elongated at a high confidence level. DAOFIND, 
the task which identifies potential objects in the images, 
only identifies point sources based on the presence of a 
distinct peak. Thus, automated photometry will be bi- 
ased against the detection of very close binaries. This 
limitation can be overcome for known or suspected bi- 
naries by manually adding a second point source in ap- 
proximately the correct location and letting ALLSTAR 
recenter it to optimize the fit; if it does not produce a 
statistically significant fit, ALLSTAR then discards it. 
We have done this for USco-109. 

In Figure 2, we present plots for USco-109 and the 
next-brightest object in the same images, a likely back- 
ground field star (Section 3.2). The first three columns 
show USco-109 and the residuals from fitting with one 
and then two point sources, and the last two columns 
show the field star and its single-source fit. The maxi- 
mum and minimum pixel values are also given to allow 
quantitative comparison of the residuals to the original 
images. The common position angle of the residuals in 
the single-source fit in all three filters seems to imply 
that this elongation is a real effect, and not simply noise. 
USco-109 appears elongated in the same direction rela- 
tive to both its neighbor and USco-67; this extension is 
therefore unlikely to be an artifact due to excess jitter, 
which would affect the field star as well. USco-109 is 
better fit with two point sources in all three filters, and 
the fit reduces the residuals by factors of 2-8. Moreover, 
these fits independently find similar positions and flux 
ratios in all three filters, which further supports its classi- 
fication as a binary system. Since the similar two-source 
residuals in each filter imply some remaining uncertainty 
in the fit and the separation is only ~1.3 pixels, foUowup 
observation of this system to confirm its multiplicity and 
properties should be a priority. However, we will pro- 

The V flux for USco-112 appears to be anomalously bright. 
Based on its i' and z' magnitudes, it appears to be of similar bright- 
ness and color to USco-55 B. However, it is 0.88 ± 0.06 magnitudes 
brighter than USco-55 B in V . Exposures in both dither positions 
give consistent fluxes, so it is probably not the result of a cosmic 
ray hit on one exposure. We suggest that this is caused either by 
a transient optical brightening of the system (e.g. a flare or in- 
creased accretion) or by contamination from a spatially unresolved 
WiiPbackground object. 
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TABLE 1 

Very Low Mass Objects in Upper Scorpius 



Name 


Date'' 




■/V-i 




SpT'^ 


M (Mq)'' 


USco-55 A 


49.9 


18.73±0.04 


15.92±0.012 


14.81±0.09 


M5.5 


0.10±0.03 


USco-55 B 


49.9 


18.92±0.015 


16.19±0.08 


15.12±0.09 


{M6.0) 


0.07±0.02 


USco-66 A 


62.6 


18.92±0.021 


16.36±0.016 


15.41±0.012 


M6.0 


0.07±0.02 


USco-66 B 


62.6 


18.94±0.022 


16.30±0.023 


15.29±0.008 


(M6.0) 


0.07±0.02 


USco-67 


49.0 


18.47±0.002 


15.52±0.001 


14.32±0.001 


M5.5 


0.10±0.03 


USco-75 


59.7 


18.71±0.002 


15.80±0.001 


14.61±0.001 


M6.0 


0.07±0.02 


USco-100 


67.0 


19.09±0.002 


15.99±0.001 


14.73±0.001 


M7.0 


0.05±0.01 


USco-109 A 


64.6 


20.20±0.03 


17.13±0.04 


15.87±0.016 


M6.0 


0.07±0.02 


USco-109 B 


64.6 


21.21±0.076 


18.07±0.08 


16.92±0.04 


{M7.5) 


0.04±0.01 


USco-112 


71.6 


18.04±0.001° 


16.15±0.001 


15.17±0.001 


M5.5 


0.10±0.03 


USco-128 


70.6 


21.29±0.007 


17.79±0.002 


16.34±0.002 


M7.0 


0.05±0.01 


USco-130 


87.8 


21.39±0.008 


17.80±0.002 


16.30±0.002 


M7.5'' 


0.04±0.01 


USco-131 


63.6 


21.98±0.011 


18.22±0.002 


16.67±0.002 


M6.5'' 


0.06±0.01 


USco-132 


88.8 


22.11±0.013 


18.31±0.002 


16.59±0.002 


M7.0 


0.05±0.01 


USco-137 


88.9 


22.92±0.023 


19.27±0.004 


17.75±0.003 


M7.0f 


0.05±0.01 



''Observation Date: JD minus 2452800. 

^Uncertainties are statistical; systematic uncertainties due to aperture corrections and conver- 
sion to standard systems are ~0.03 magnitudes. 

'^Spectral types in parentheses are inferred from photometry presented here; others are from 
Ardilla et ah(2000) (Section 4.1) 

^iVIasses are determined from the models of Baraffe et al. (1998). 

^The V magnitude for USco-112 is anomalously bright. (Section 2) 

^Spectral type is based only on photometric R ~ Iq colors. 



ceed under the assumption it is a binary system in our 
subsequent analysis. 

We summarize the properties of these three systems 
in Table 2. The uncertainties in separation and position 
angle are determined from the standard deviations in the 
locations as reported by ALLSTAR for the three filters. 

3.2. VLMOs and Background Stars 

In Figure 3, we present an i' versus i' — z' color magni- 
tude diagram of the 12 primary targets (filled circles) and 
all other objects detected at the 5a level in both filters. 
Objects within 5" (725 AU), which are statistically more 
likely to be physically associated, are shown as open cir- 
cles while other objects are shown as crosses. Also shown 
are t he average main seq uence at the distance of Upper 
Sco l)Hawlev et all 1200^ and a 5 Myr iso chrone based 
on the evolutionary models of lBaraffe et a l. (1998). The 
location of the isochrone is determined by converting the 
predicted Ic and J magnitudes to i' and z' magnitudes, 
using Ic — i' colors derived using the methods described 
in Section 2 and z' — J colors found from the SDSS 
fiel d main sequence for older, more massive M dwarfs 
by 'Hawle v et all lj2002(l . Although the Ic - i' trans- 
formations should be accurate (Section 2), the z' — J 
transformations may be more sensitive to surface gravity 
differences due to the larger difference in central wave- 
lengths. Consequently, the i' — z' color of the 5 Myr 
isochrone is somewhat uncertain. 

The three very close companions to USco-55, USco-66, 
and USco-109 are located above the SDSS main sequence 
and well above the background population, so they are 
most likely association members. Based on their close 
proximity and the low density of association members, 
we conclude that these are physically associated com- 
panions. All other objects fall well below the empirical 
main sequence and are most likely background stars. The 
marginally resolved companion to USco-109 resides near 



the empirical (but somewhat uncertain) main sequence, 
but given the uncertainty in its i' — z' color, its location 
is consistent with the broad luminosity distribution of 
low-mass Upper Sco members that we discuss in Section 
5. 

3.3. Sensitivity Limits 

We determined detection limits as a function of dis- 
tance from the primar y stars via a Mo nte Ca rlo simula- 
tion similar to that of iMetc hcv et al] l)2003D . We used 
the IRAF task DAOPHOT/ADDSTAR and our average 
PSF to add artificial stars at a range of radial separations 
and magnitudes to the fields of USco-67, USco-128, and 
USco-137, which represent the maximum, median, and 
minimum brightness sources in our sample. We then at- 
tempted to identify them with our ALLSTAR photome- 
try pipeline. In Figure 4, we show the detection limits for 
the USco-67 and USco-128 fields in V and z', as a func- 
tion of separation, at which we can detect > 10%, 50%, 
and 90% of the companions. At small separations (<20 
AU) , the 50% detection thresholds roughly scale with the 
brightness of the primary; they are similar for all objects 
in terms of Am. The detection limits converge to con- 
stant values at large separations where the background 
dominates the noise: z'=23.40±0.05, i' =25.08±0.05, 
and 1/=26.01±0.05 (uncertainties are based on the stan- 
dard deviation between the three fields). The simulations 
demonstrate that we potentially could identify bright 
pairs as close at 1 pixel (0.027"; 4 AU) but none closer. 
We also show the locations (in Am and separation) of 
the three companions to USco-55, 66, and 109. These 
results indicate that the probability of detecting a sys- 
tem with the separation and flux ratio of USco-109 with 
automated photometry is between 10% and 50% for the 
brightest objects and becomes negligible for the fainter 
targets in our sample. However, as discussed in Section 
2, these limits are conservative at small separations since 
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TABLE 2 

Binary Parameters 



Name 


Projected 


Projected 


Position 


ASpT 


g 




Separation (mas) 


Separation(AU) 


Angle(deg) 






USco-55 AB 


121.6±0.6 


17.63±0.09 


307.7±0.4 


O.SitO.l 


0.86±0.04 


USco-66 AB 


70.3±0.5 


10.19±0.07 


31.7±0.2 


O.OitO.l 


0.94±0.03 


USco-109 AB 


34±2 


4.9±0.3 


302±3 


1.5±0.2 


0.59±0.04 



ALLSTAR often reports only a single point source when 
the PSF is still noticeably elongated and manual fitting 
can distinguish two point sources. The position angle of 
the binary can also play a role in the detection proba- 
bility; as can be seen in Figure 1, the blended PSF of 
USco-109 is elongated along the axis where the PSF is 
narrowest, providing the most significant possible detec- 
tion. Artificial star tests show that if the companion was 
near one of the extended lobes of the PSF, it would have 
been much more difficult to identify. 

3.4. Uncertainties in Binary Properties 

A similar Monte Carlo routine was used to test the un- 
certainties in the measurements of USco-109. We used 
ADDSTAR and the average PSF to construct 100 sim- 
ulated images, given the positions and brightnesses re- 
ported for the primary and the secondary for the real 
images, and then used ALLSTAR to perform PSF-fitting 
photometry on these simulated images. The standard de- 
viation in separation ('--^l mas) from the simulated images 
is consistent with that calculated from the standard de- 
viation in separation between the three filters (~2 mas), 
and the standard deviation in the flux ratio Am (0.09 
magnitudes in V, 0.06 in i' and z') is consistent with 
that determined from the magnitudes reported by ALL- 
STAR for each object. 

4. ANALYSIS 
4.1. Inferred Properties 

In Table 1, we give the inferred spectral types and 
masses for all of the VLMOs in our sample. Spec- 
tral types for single VLMOs and for binary pr imaries 
in our sample are taken from lArdila et all lj200 0fl. Spec- 
tral types for nine objects were determined from low- 
resolution spectra; the other three (USco-130, 131, and 
137) have only photometric spectral types as determined 
from R—Ic colors. Ardila et al. found that although this 
region is young, the extinction is generally low [Ay < 1), 
so the photometric spectral types should be reasonably 
accurate. Sp ectroscopic obse r vation s of all targets except 
USco-137 bv lMohantv et all 1)2005(1 appear to be consis- 
tent with these assigned spectral types. 

The masses for this sample are estimated from 
the 5-Myr m ass-magnitude-temperature relations of 
iBaraffc et al. (1998) and the temperature-SpT relations 
o fiLuhman et al. 1 2003) . and range from 0.04 to 0.10 Mq. 
Large systematic errors may be present in these and all 
pre-main sequence models (e.g. Baraffe et al. 2002; Close 
et al. 2005; Reiners et al. 2005), so they are best used 
for relative comparison only. The models of Baraffe et 
al. predict that the mass ratio (q = rris/mp) and ASpT 
for the three binaries are a function of the primary-to- 



secondary flux ratios Am with only a minor mass de- 
pendence; we report these quantities as determined from 
the flux ratio Ai' in Table 1. The uncertainties reflect 
the uncertainties in the flux ratios, but do not include 
any systematic uncertainties from the models. The val- 
ues determined from AV^ and Az' are consistent with 
these results, but the values for Az' are least dependent 
on primary mass, so they are most reliable. 

4.2. The Binary Frequency at the Stellar /Substellar 
Boundary 

For our 12 M5.5-M7.5 targets, we have detected three 
binary VLMO systems with projected separations of 4.9- 
17.6 AU, giving an observed binary fraction of 25±14%. 
However, as can be seen in Figure 4, the detection of faint 
companions is difficult at separations comparable to the 
PSF width (58 mas in V). Consequently, the total binary 
fraction may be higher. We can estimate the number of 
undetected companions via a method based on that of 
iClose e t al. (2003) where we convolve an assumed mass 
ratio distribu t ion w ith the mass-magnitude relations of 
IBaraffe et alJ l)1998j) and the detection thresholds Am as 
calculated in section 4.3. Similar results are found for 
all three filters; the following results are those for the i' 
band. 

Observations of field M-dwarfs have found a flat mass 
ratio distribution from 1.0 to 0.1 l)Fischer fc Marcvl 1992(1 

for all companions at separations a > 4 AU. However, 
studies of field VLMOs have found only BD companions 
with mass ratios q > 0.6 (Am < 1). If we adopt a mass 
ratio distribution that is uniform between 0.6 and 1.0 
and zero for q < 0.6, then the detection thresholds in 
Figure 4 can be used to predict the number of unidenti- 
fied companions. The 50% detection limit spans the full 
range of possible mass ratios at separations >10 AU, and 
then approaches Am = at a separation of ^4 AU. The 
average threshold of Am = 0.5 in the separation range 
4 < a < 10 AU corresponds to a detection threshold of 
q = 0.8. This implies that all companions outside 10 AU 
and ^1/2 of all companions inside 10 AU were detected. 
Since we detected one companion within 10 AU, we pre- 
dict the existence of one additional companion, which 
yields a completeness-corrected binary fraction for sep- 
arations >4 AU of 33 ± 17%. If we assume a flat mass 
ratio distribution between 0.1 and 1.0, as for field stars, 
then similar analysis finds 1/2 of an undetected compan- 
ion for separations 10 < a < 15 AU and 2 companions for 
a < 10 AU, for a complete binary fraction of 46 ± 25%. 
However, the paucity of observed unequal-mass compan- 
ions strongly implies that the lower value is more appro- 
priate, and even it might be an overestimate since our 
Monte Carlo tests produced conservative limits at small 
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separations and USco-109 lies below our nominal detec- 
tion curve. 

We note that the intrinsically higher luminosity of bi- 
naries makes them easier to identify in membership sur- 
veys, which can bias the binary frequency to larger val- 
ues. However, this should not significantly bias our re- 
sults since all targets except USco-137 are substantially 
brighter than the detection limit of the original survey 
by Ardila et al. 

4.3. Limits on Planetary- Mass Companions 

The high dynamic range observations of young stars 
and brown dwarfs presented here have the potential to 
directly image wide planetary-mass companions. In Fig- 
ure 4, we indicate the predicted zero-extinction bright- 
ness of some representative masses of planetary compan- 
ions based on the models of Burrows et al. (1997); since 
most of our targets have Ay < 1, extinction correc- 
tions at i' and z' should be small (< 0.5 magnitudes). 
The extremely red colors predicted for planetary-mass 
objects imply that z' observations provide the strictest 
limits on planetary companions. Based on the lack of 
detections, we conclude that there are no planetary com- 
panions with mass >5 Mj at projected separations larger 
than 280 mas (40 AU) or mass >10 Mj at projected sep- 
arations larger than 140 mas (20 AU) among our sam- 
ple. For comparison, these limits would have allowed for 
marginal detection of the candidate planetary compan- 
ion to 2MASSWJ 1207334-393254, a substellar member 
of the TW Hya association, which has a projected sepa- 
ration of ~54 A U, a flux ratio of A^:'^ 7, and a predicted 
mass of -5 Mj l|Chauvin et al.ll2004|) . 

5. DISCUSSION 

Multiplicity surveys of VLMOs have suggested a 
fairly uniform set of binary proper ties in th e field 
llClose all I2003t iBouv et all 120031: lBurgassc r~etaLl 
12003': 'Sieder ct al I2005D and in the Pleiades 
l|Mar tin et al. 2003). No companions with wide 
separations (>20 AU) or with unequal mass ratios 
{q < 0.7) were found, despite sufficient sensitivity for 
their detection. The binary fractions observed were also 
significantly lower than the binary fractions of 57% for 
field G-dwarfs (Duauennov & Mavor 199D and 35-43% 
for fi eld M-dwarfs ( Reid & Gizis 199lil iFischer fc Marcvl 
Il992j) . Specifically, Close et al. found a binary fraction 
of 15 ± 7% for separations of >3 AU for field late-M 
and L dwarfs, Bouy et al. found a binary fraction of 
18.8 ± 3.7% at separations of >1 AU for field L dwarfs, 
Burgasser et al. found a binary fraction of 9^g^% at 
separations of >1 AU for field T dwarfs, Siegler et 
al. found a binary fraction of 9^3% for field M6-M7.5 
dwarfs at separations of >3 AU, and Martin et al. found 
a binary fraction of 15^5^% for separations of >7 AU 
for Pleiades BDs. 

The study presented here is the first survey of young 
BDs with the resolution and sensitivity to identify bina- 
ries as close as 4 AU, allowing a more direct compari- 
son with field surveys. Our method of correcting for de- 
tection biases is based on that of Close ct al. (200^, so 
our completeness-corrected binary fraction can be com- 
pared most directly to that study. Our binary fraction, 
33 ± 17% for separations >4 AU, is higher than, but 



statistically consistent with their completeness-corrected 
binary fraction for field VLMOs (15 ± 7%). However, 
the field surveys are sensitive to companions with sepa- 
rations as small as ~ 1 AU and found many companions 
within 4 AU that we could not have detected because of 
the larger distance to our targets. Our binary fraction for 
companions with a > 4 AU only sets a lower limit on the 
total fraction of VLMO binaries. Higher-resolution ob- 
servations or long-term spectroscopic monitoring of these 
young VLMOs is needed to determine if they too have 
an abundance of even closer companions. 

All three Upper Sco VLMO binaries have flux ratios of 
Ar7i<1.0, corresponding to mass ratios of > 0.6, which 
also agrees with with the other surveys described above, 
which found no BD binaries with q < 0.6. The implica- 
tion is that the mass ratio distribution is not flat at low 
masses, but instead is biased toward binaries of nearly- 
equal mass. This indicates that a completeness correc- 
tion which assumes a flat mass ratio for q^O.l would 
strongly overestimate the VLMO binary frequency. Like 
the field and open cluster surveys, we found no wide bi- 
nary companions. Our results are also consistent with 
a surve y of the y oung (^3 Myr) star-forming region IC 
348 bv iLuhmanI l|20 04b ) . which found a Icr upper limit 
on the binary fraction of ~7% for separations >20 AU. 
Wide (> 20 AU) companions are therefore much less 
common around VLMOs than around G and early M 
dwarfs, though some candidate wide companions have 
been identified in Chamaeleon I l lNeuhauser et al.]l2002l: 
lLuhmanll2004ah and TW Hydra l|Chauvin et alJl2004D . 

The distinct binary properties of field BDs, relative to 
those of stars, indicate that these properties are mass 
depend ent. However, the form of this dependence is not 
known: iKroupa et alJ (|2003) and lClose et~an l|2003|) in- 
terpret current data as a break near the stellar/substellar 
boundary while Luhman (2004b) argues for a smooth 
mass dependence. We investigate the mass dependence 
at a young age by combining our results with a binary 
ce nsus of more massive T Tauri members of Upper Sco 
bv IKohler et alJ lj2000|) . This survey targetted 118 X- 
ray selected pre-main-sequence stars ranging from G5 to 
M5, corresponding to masses of 2.0 — O.13M0 and found 
a binary fraction of 35.2 ± 6.3% for companions with 
separations between 6" (870 AU) and 0.13" (19 AU). 
Dividing the sample into spectral type bins reveals no 
spectral type dependence of the binary frequency, though 
the subsample sizes are not sufficient to completely rule 
out a gradual decline. Specifically, the survey found a 
binary frequency of 26/50 = 52 ± 10% for Upper Sco 
members of spectral type M0-M5 (0.7 - 0.13Mo) for 
separations between 19 and 870 AU, including several 
very wide companions to primaries in the M4-M5 range. 
Kohler et al. did not confirm association membership 
of the companions with colors or proper motions, but 
they estimate from the low background object density 
(6.64 ± 0.45 deg~^ in K, comparable to our background 
object density of 7.8 ± 0.8 deg''^ in i' and z') that no 
more than 15% of their detected binaries were spurious 
claims. They also found no decline in mean or max- 
imum separation and no preference toward equal-mass 
companions with decreasing mass, though the mass ra- 
tio result could not be tested for the lowest-mass targets 
since unequal-mass companions would have fallen below 
the detection limit. Our measurements in the mass range 
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0.04 — 0.1 OM0 suggest a significant decrease in the num- 
ber of wide binaries. Since none of the newly-discovered 
binaries presented here would have been found in the 
Kohler survey, there may be a discontinuity in the sepa- 
ration distribution at a mass of '^O.IOM^q. 

We also note that resolving the newly-discovered bi- 
naries does not reduce the large scatter in the color- 
magnitude diagram at low mas ses that has b een ob- 
served in previous surveys ijArd ila et al] 120001) and is 
evident for the app arently single VLMOs in Figure 1. 
iPreibisch fc Ziniiecker ( 199S.) argue for a consistent age 
of 5 Myr based on the narrow locus seen on color- 
magnitude diagrams of more massive Upper Sco mem- 
bers, but the p oor fit between our data and any sin- 
gle isochrone of iBaraffe et al.l (jlOOS^) indicates that the 
lowest-mass members of Upper Sco could have a signif- 
icant range in ages or distances. The uncertainty intro- 
duced by transforming of the 5 Myr isochrone into the 
SDSS magnitude system makes it difficult to interpret 
other possibilities, such as the existence of additional 
unresolved binaries or preferential selection of older or 
younger systems in the original survey by Ardila et al. 

In summary, the binary properties of the young VL- 
MOs observed here are the first convincing evidence that 
the small separations and equal mass ratios of VLMO bi- 
naries are established at young ages rather than as the re- 
sult of subsequent dynamical evolution after the T Tauri 
phase of evolution. The uncertainty in our total binary 



fraction is too large to directly test the suggestion that 
the binary fraction declines with declining mass, but if 
the total binary fraction of the field population is truly 
primordial, as the small separations and equal masses 
appear to be, then this provides powerful evidence in 
support of this assertion. However, complementary spec- 
troscopic surveys are needed to determine if the apparent 
deficit of substellar binaries is a true absence or is due to 
a shift to smaller separations than can be detected with 
imaging alone. Finally, the recent discovery of candi- 
date wide binary systems in the low-densi ty star-forming 

2002; LuhmaJ 



region s Chamaeleon 1 ^Ngyhgjjgg^^l 



I2004a|) and TW Hydra ijChauvin et al 



20041 indicate 



that there could be differences in the VLMO binary for- 
mation process and initial conditions in OB associations 
compared to low-density T associations. 
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USco 55 USco 66 USco 109 USco 67 (single) 




pixels 



Fig. 1. — Contour plots of four targets: USco-55, USco-66, USco-109, and USco-67, respectively. Units are in pixels, and the projected 
physical scale at the distance of Upper Sco is shown in the upper right panel. Contours arc drawn at 95% through 5% of the mEiximum 
pixel value, in increments of 10%. The field of view in each image is 432 mas, or "^60 AU at 145 pc. 
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USco 109 Single Fit Double Fit Field Star Single Fit 




pixels 



Fig. 2. Contour plots of USco-109 and a background star in the same field for all three filters. The first three eoluiniis show USco-109 
and the residuals from fitting with one and then two point sources, and the last two columns show the background star and the residuals 
from a single-source fit. For residuals, contours arc drawn at the 90%, 50%, and 10% levels of maxinmm (solid lines) and minimum (dashed 
lines). The maximum and minimum pixel values arc given to allow comparison of the residuals to the original images. The pixel values in 
the last column, where the sky background contours fill each panel and obscure the text, are (58,-133), (29,-27), and (10,-9), respectively. 
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0.5 1 1.5 

i' — z' (mag) 

Fig. 3. — An i' vs i' — z' color magnitude diagram. Primary Upper Sco targets are shown as filled circles, neighbors within 5" arc shown 
as open circles, and widely-separated objects are shown as crosses. The SDSS field main sequence (solid line), a 5 Myr isochrone (dashed 
line), and the detection limits of the survey (dotted line) are also shown. Candidate binary VLMO pairs are connected with dotted lines 
and labeled. The error bars for binary components are associated with each point, and the error bars for single objects are shown in the 
upper left corner. 
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Fig. 4. — Five-sigma detection frequencies (10%: solid, 50%: dashed, and 90%: dash-dotted) as a function of separation in each filter 
for the VLMOs of maximum and median brightness, USco-67 and USco-128. Corresponding brightnesses of potential planetary-mass 
companions are shown on the right for the z' plot. The brightness of the primary object is denoted with a star to allow conversion to Am 
values, and the vertical dotted line indicates a separation of 20 AU at the distance of Upper Sco (145 pc). The filled circles mark the 
separation and Am values for the three binary systems. 



